Gene functions of the Ambystoma altamirani skin microbiome vary across space and time but potential antifungal genes are widespread and prevalent.

Amphibian skin microbiomes can play a critical role in host survival against emerging diseases by protecting their host against pathogens. While a plethora of biotic and abiotic factors have been shown to influence the taxonomic diversity of amphibian skin microbiomes it remains unclear whether functional genomic diversity varies in response to temporal and environmental factors. Here we applied a metagenomic approach to evaluate whether seasonality, distinct elevations/sites, and pathogen presence influenced the functional genomic diversity of the A. altamirani skin microbiome. We obtained a gene catalogue of 92 107 nonredundant annotated genes and a set of 50 unique metagenome assembled genomes (MAGs). Our analysis showed that genes linked to general and potential antifungal traits significantly differed across seasons and sampling locations at different elevations. Moreover, we found that the functional genomic diversity of A. altamirani skin microbiome differed between B. dendrobatidis infected and not infected axolotls only during winter, suggesting an interaction between seasonality and pathogen infection. In addition, we identified the presence of genes and biosynthetic gene clusters (BGCs) linked to potential antifungal functions such as biofilm formation, quorum sensing, secretion systems, secondary metabolite biosynthesis, and chitin degradation. Interestingly genes linked to these potential antifungal traits were mainly identified in Burkholderiales and Chitinophagales MAGs. Overall, our results identified functional traits linked to potential antifungal functions in the A. altamirani skin microbiome regardless of variation in the functional diversity across seasons, elevations/sites, and pathogen presence. Our findings suggest that potential antifungal traits found in Burkholderiales and Chitinophagales taxa could be related to the capacity of A. altamirani to survive in the presence of Bd, although further experimental analyses are required to test this hypothesis.

The skin microbiota of the axolotl Ambystoma altamirani is highly influenced by metamorphosis and seasonality but not by pathogen infection.

BACKGROUND: Microbiomes have been increasingly recognized as major contributors to host health and survival. In amphibians, bacterial members of the skin microbiota protect their hosts by inhibiting the growth of the fungal pathogen Batrachochytrium dendrobatidis (Bd). Even though several studies describe the influence of biotic and abiotic factors over the skin microbiota, it remains unclear how these symbiotic bacterial communities vary across time and development. This is particularly relevant for species that undergo metamorphosis as it has been shown that host physiology and ecology drastically influence diversity of the skin microbiome. RESULTS: We found that the skin bacterial communities of the axolotl A. altamirani are largely influenced by the metamorphic status of the host and by seasonal variation of abiotic factors such as temperature, pH, dissolved oxygen and conductivity. Despite high Bd prevalence in these samples, the bacterial diversity of the skin microbiota did not differ between infected and non-infected axolotls, although relative abundance of particular bacteria were correlated with Bd infection intensity. CONCLUSIONS: Our work shows that metamorphosis is a crucial process that shapes skin bacterial communities and that axolotls under different developmental stages respond differently to environmental seasonal variations. Moreover, this study greatly contributes to a better understanding of the factors that shape amphibian skin microbiota, especially in a largely underexplored group like axolotls (Mexican Ambystoma species).